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In this paper we consider the possible existence of unstable axisymmetric modes in Kerr space 
times, resulting from exponentially growing solutions of the Teukolsky equation. We describe a trans- 
formation that casts the radial equation that results upon separation of variables in the Teukolsky 
equation, in the form of a Schrodinger equation, and combine the properties of the solutions of this 
equations with some recent results on the asymptotic behaviour of spin weighted spheroidal harmon- 
ics to prove the existence of an infinite family of unstable modes. Thus we prove that the stationary 
region beyond a Kerr black hole inner horizon is unstable under gravitational linear perturbations. 
We also prove that Kerr space-time with angular momentum larger than its square mass, which has 
a naked singularity, is unstable. 

PACS numbers: 04.50.+h,04.20.-q,04.70.-s, 04.30.-w 
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I. INTRODUCTION 



In 1916, shortly after Einstein published his vacuum field equations, Schwarzschild found a solution that represents 
what came to be known as a static black hole. It took almost fifty years to find a metric representing a rotating black 
hole. This is Kerr's solution [![, which, in Boyer-Lindquist coordinates reads 



ds 2 = - 



(A - a 2 sin 2 



-dt - 2a sin 



(r 2 + a 2 - A) 



dtdcf) - 



(r 2 + a 2 ) 2 - Aa 2 sin 2 i 



sin- ti(nr + —dr 2 + Ed6> , (1) 



where E = r + a cos 9 and A = r — 2 Air + a . Kerr's metric has two integration constants: the mass M, and 
the angular momentum per unit mass a. We will only consider the case M > 0, and we will take a > without 
loss of generality, since for a < we can always change coordinates <f> — > — 4>, under which a — > —a. Setting a = 
(no rotation) gives Schwarzschild's solution. The E = ring singularity at r = 0,9 — ir/2 is hidden behind a black 
hole horizon as long as a < M. The inner and outer horizons are located at the zeroes of A: = M — V M 2 — a 2 
and r = M + V M 2 — a 2 , which are just coordinate singularities in ([1]). Kerr's space-time can be extended through 
these horizons, and new regions isometric to / : r > r a , II : ri < r < r a and III : r < Ti arise ad infinitum, giving 
rise to the Penrose diagram displayed in Figure 1. According to the Carter- Robinson theorem and further results 
by Hawking and Wald, (see, e.g., [2| and references therein) if (A4,g a b) is an asymptotically-flat stationary vacuum 
black hole that is non-singular on and outside its event horizon, then it must be a < a < M member of the 
two-parameter Kerr family. In other words, the metric (JTJ) with a < M is the unique rotating black hole solution of 
Eintein's equations. All these features, together with its relevance to model astrophysical processes, signal Kerr's 
metric as one of the most important among the known exact solutions of Einstein's equations. 



Regions I and III of a Kerr black hole are stationary, and finding out whether or not they are stable under 
gravitational perturbations is a fundamental issue. The stability of region I was established after a series of papers 
following Teukolsky's discovery of the perturbed fields admitting separable evolution equations Q , and ending with 
Whiting's proof that no solution of Teukolsky equations, supported in region I and properly decaying at infinity and 
r , can grow exponentially in time [J]. According to these results, small disturbances constrained to region I will 
oscillate around ([T]), in the same way a mechanical system oscillates around a local minimum configuration of its 
potential energy. Region III has a number of undesirable properties: it admits closed time like curves, has the ring 
singularity, and lies beyond the Cauchy horizon, so that predictability is lost for earlier data. It has been argued Q 
that a "mass inflation" mechanism would prevent the formation of the Cauchy horizon, that would then be replaced 
by a null singularity acting as a boundary of the spacetime and cutting off this region. A fundamental issue to address 
is then the analysis of the linear stability of this region under gravitational perturbations. In this paper we show that 
region III is unstable by proving that there are solutions of Teukolsky equations, restricted to region III and satisfying 
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FIG. 1: Penrose diagram for the maximal analytic extension of Kerr's space-time. Regions labeled I and I' are isometric, and 
so are II and II', and III , III' 

appropriate boundary conditions, that grow exponentially in time. This instability of Kerr space-time remained 
uncovered since Kerr's metric discovery forty five years ago and makes the analytic extension of Kerr's solution less 
relevant. One is immediately led to ask where is region III driven to by a gravitational perturbation that excites 
the unstable modes. Unfortunately, answering this question goes beyond the linear perturbation analysis in this paper. 

We also prove in this paper that Q]) with a > M is unstable. For this "super extremal" case, the ring singularity at 
r = 0, 8 = 7r/2 is naked, i.e., not protected by a black hole horizon, and can therefore causally connect with the rest 
of space-time. The super extremal case was considered recently in Q , where numerical evidence for the existence of 
unstable modes was presented. Because of the nature of the numerical procedures used, the evidence was restricted 
to a range of values of a/M. Clearly, a more general result, establishing in a rigorous way the unstable nature of the 
Kerr space time in the case a > M is of interest, as it suggests that such space times cannot be the end point of the 
collapse of a self gravitating system, a result of relevance in relation to the Cosmic Censorship Conjecture. For other 
nakedly singular space times such as the negative mass Schwarzschild solution and super extreme Reissner-Nordstrom 
spacetimes, the linear instability was established after 0, [1] and [1, Q respectively, where explicit expressions for 
unstable modes can be found. It should be stressed, however, that a > M Kerr space-times are a subset of a larger 
family of axially symmetric solutions with a naked singularity, not singled out by the uniqueness theorems that make 
the black hole case so relevant. 

This paper is organized as follows: Section [TT] reviews Teukolsky radial and angular equations, focusing on axisym- 
metric unstable modes. Although all we need from the angular equation is the behaviour of eigenvalues for small 
and large purely imaginary frequency, available in the recent literature [l2| , a technique to obtain approximations for 
eigenvalues and eigenfunctions is developed for completeness, and tested against the results in [12fl . The reformulation 
of the radial Teukolsky equation that follows is crucial in this work. The radial equation is cast in Schrodinger like 
form, and its potential analyzed. This formulation is used in Section IIIII to prove that there exist infinitely many 
unstable modes for the super-extremal Kerr spacetime. Finally, in Section ITVl we adapt this proof to show that region 
III of a Kerr black hole is also unstable. 
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II. TEUKOLSKY EQUATIONS 

The systematic analysis of the linearized perturbations of Kerr space-time was greatly facilitated by Tcukolsky's 
discovery [3[ that, in the Newman-Penrose formalism, one can write second order partial differential evolution equations 
for two linearized null tetrad components <I> S of the Weyl tensor, with s = ±2 the spin weight of these components. 
Furthermore these equations are separable. Introducing 

<S>s = R w , m , s {r)S™ s (d)exp{im<f>)exp(-iLot), (2) 

the linearized equations reduce to a coupled system for S and R, 

1 d ( . AS\ / 9 9 9 „ „ (m + scos9) 2 _ 9 \ _ 

^ sin 0—\ + [a 2 u 2 cos 2 9 - 2atos cos 9-- 5 '- + E - s 2 ) S = (3) 

smOdO V dO \ sin 2 6 1 



d 2 R dA dR ( K 2 - 2is(r — M)K , 2 



A— r + (s + 1)- - + < j h 4irws — [E — 2amu + a 2 tu 2 - s(s + l)}> R = 0, (4) 

dr z dr dr [ A J 

where A = r 2 — 2Mr + a 2 , and K = (r 2 + a 2 )to — am, with M the mass, and a the rotation parameter, related 
to the angular momentum J by J = aM. Tcukolsky's equations also apply to scalar (s — 0) and electromagnetic 
(s = ±1) perturbations on the Kerr background. Moreover, any solution ip of the s = —2 Tcukolsky equations is 
a "Debye potential" from which a metric perturbation h a b that solves the linearized Einstein equations around the 
Kerr b ackg round can be constructed, h a b being obtained by applying a second order linear differential operator to 



tjj (see [13j). Teukolsky's equations have been applied in a wide range of problems related to the behavior of Kerr 
black holes under gravitational and other perturbations, including black hole collisions in the so called "close limit" 
[lij . A very important question addressed to from the early applications is that of the stability of region I of a Kerr 
black hole under gravitational perturbations, which was established through a series of works starting from 0] and 
ending with 0]. In what follows we will prove that region III of an a < M Kerr spacetime, as well as the nakedly 
singular a > M Kerr spacetimes, are unstable. This is done by showing that for these cases there exist solutions 
of the Teukolsky equations with purely imaginary frequency uj and satisfying appropriate boundary conditions. Our 
proof is based on an analysis of the behavior of the lowest eigenvalues E of ([3]) and (|4]) for large purely imaginary 
frequency to. Note that the coupling in the system ©-((H) is rather subtle, and comes from this common eigenvalue 
E. Suppose s and m are fixed, then given u>, E has to be chosen so that S is regular on the sphere (such a solution 
is called a spin weighted spheroidal harmonic, SWSH), and such that R satisfies the required boundary conditions. 
The first condition gives a countable set of possible values Ei(auj), whereas the second one fixes the possible w values. 
Each one of these conditions can be cast in the form of a continued fraction equation [Icl ] that arises when solving 
the three term recursion relation on the coefficients of series solution for S and R in ^ and (UJ). 

A. Spin weighted spheroidal harmonics 



For complex u> near u> — 0, a Taylor expansion for Ei(auj) up to order (auj) 6 is available in the literature ( 111. 



[l2| and references therein) . Asymptotic expansions for to — > oo and for to — > +ioo can be found in [l2l [TEl . flfl 
Particularly useful to our purposes is the fact that in the axial case m = 0, and for gravitational perturbations s = ±2, 
the E eigenvalue for large purely imaginary to behaves as fl2j 

B/HU=i* = (2e-3)k + O(k ), asfc^oo (5) 

In the above formula it was assumed that that the I labeling is chosen such that the (I, m, s) SWSH reduces to the 
corresponding spin weighted spherical harmonic in the ato = limit [12j , which implies 

Ee(aLo)\ au=0 =£(£+l), (6) 

the relevant modes being £ = 2, 3, 4, ... 

For the purposes of this paper we set m — 0,au> = ik, s = —2, and, after changing variables to x = cos 9 in ©, we 
find that the angular Teukolsky equation is 



d 
dx 



(l-x 2 )f 
dx 



k 2 x 2 + )S + i AkxS = -ES. (7) 

1 — x z 



If k = 0, reduces to an associated Legendre equation, whose regular solutions are the associated Legendre polyno- 

(2) 

mials Pi {x),£ — 2,3, for which ^ holds. For real nonzero k, |7]) is a complex equation, but still the allowed E 
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FIG. 2: The left panel shows the lowest (corresponding to I — 2) E roots obtained by solving eq. ((9J for thirty different values 
of k £ [0, 15], and a least square linear fit for the ten points in the interval 10 < k < 15, which gave E = 0.993fc+constant, in 
excellent agreement with the expected asymptotic expansion Eq. ©. The right panel shows forty values obtained in this way 
for < k < 4 (dots) and the low frequency approximation in [Tl|] [12| ] up to order k 4 (solid line). 



values are real (see e.g., equation (2.3) in [12]). The symmetries of ([7|) for real E imply that if S(x) is a solution then 
so are S(—x) and S±(x) := S(x) ± S(— x) (an overline denotes complex conjugation), so we can restrict our search to 
solutions with real and imaginary parts of opposite parity and, multiplying by i if necessary, we can always assume 
that S(x) has even real and odd imaginary parts [22j: 

S(x)^^x) (8) 

Equation |(7J) has regular singular points at x — ±1. There is a regular Frobenius solution around x = — 1 of the 
form !F(_u(x) = J2 n >i c n{E,k) (x + 1)™, and a linearly independent solution diverging as ln(a; + 1). The series in 
the regular solution converges for \x + 1| < 2 fl8|> and is finite (in fact, vanishes) at x = 1 if E = Eg(auj)\ aL j = ik for 
some I. Note that in this case !F/_x)(x) has to agree up to a multiplicative constant with the (regular) Frobenius 
series solution for © around x = 1. A SWSH is a solution of ||7J) which is regular for x € [—1,1]. Regularity 
in [—1,1] will follow from regularity at, say, x = — 1 if we impose the parity condition {§]). Note, however, that 
since S satisfies a second order ODE, the condition ((5J) is equivalent to 9(5(0)) = 5R(5'(0)) = 0. In particular, if 
5 = (^ ( _ 1} (0)) F(-i)(x) = E„,m>i c «(^> k)cm(E,k)(x + 1)™, then 3(5(0)) = 0, the solutions E{k) of the equation 

= SR(5'(0)) = (1/2) "ME, k)c m (E, k) + c n {E, k)c m {E, k)) =: V(E, k), (9) 

n,m>l 

will equal Ei(auj)\ auJ= ik for some I, and 5 will be regular for x € [—1, 1] for these E values. 

As an application of this approach, we replaced T^_i){x) in ([9]) with the regular Frobenius series at x = —1 up to 
order (x + l) 40 , and numerically solved for the lowest E (i.e., £ = 2) root in (J9j) for different values of fc. The results, 
contrasted with the low lo expansion (see [llj ] [l2j and references therein) and the asymptotic form ([5]) are shown in 
Figure 2. 



B. Teukolsky radial equation, a > M case 

Equation Q is of the form AR + QR + ZR = 0, dots denoting derivatives with respect to r, and can be put in a 
suitable form for our analysis by introducing an integrating factor L, ip := RjL and changing the radial variable to 
r* , where := 4 with / an unspecified positive definite function of r. The result is 

AL ,., {2L'A QL ALf'\ . (Ah" QL' AL'f \ , , 

~T + \F~ + ! ~F~) + (~ + f 7^ + ) = (10) 
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where ' denotes derivatives with respect to r* . The choice / = v A (recall A is strictly positive in the super extreme 
case a > M) gives an adimensional r* , 



= In 



- M + ^r 2 - 2Mr + a 2 
M 



2M\r\ 



r — > oo 
r — ► — oo 



that grows monotonically with r, and can easily be inverted in terms of elementary functions, 

Mexp(r*) „ M 2 -a 2 



M 



2Mexp(r*)' 



(11) 



(12) 



The next step is to choose L such that the tp' coefficient of (|10j) vanishes. This gives L = A 4 with (j4|) reducing 
to a Schrodinger like equation 



Hip := + VV = -Eip. 



The potential 



V := 



AL QL 



Z + E 



(13) 



(14) 



is independent of E (recall that Z is the coefficient of R in |[3])), this being the reason behind the above choice of r* 
For s — —2, m = 0, and lo — ik/a, the case we are interested in, the potential is real and reads 



V = 



r 4 + a 2 r 2 + 2a 2 Mr 
a 2 (r 2 - 2Mr + a 2 ) 



-r d + 3Mr 2 - a 2 r - a 2 M 



a(r 2 - 2Mr + a 2 ) 



r 2 - 2Mr + 15M 2 - Ua 2 



4(r 2 - 2Mr + a 2 ) 



=: rF 2 +fcVi+F 



(15) 



Since a > M in the super-extreme case, V is smooth everywhere. Also, V is bounded from below for every k > 0, its 
minimum being not continuous, as a function of k, at = 0: 



min{T/(r,A: = 0),r G M} = -7/2 
lim min{y(r, k), r G M} = -15/4 



(16) 



III. 



UNSTABLE MODES OF THE KERR NAKED SINGULARITY 



TL in (|13p is self-adjoint in the Hilbert space of square integrable functions of r* with hcrmitian product (a\f3) := 
fa/3 dr* . Given that V is smooth, bounded from below, and 



V 



Mkt 



2a 



00, 



(17) 



the spectrum of the self-adjoint operator 7i is fully discrete and has a lower bound. A careful analysis of (fTTj) . (fT5")) and 
P7p . or, more directly, the asymptotic form of indicates that the square integrable eigenfunctions of TL behave as 



■^(M) = J (l + G(M/r)) ,r->oo 
f) 1 ' (l + 0(M/r)) ,r^-oo 



(18) 



with sub-leading terms that depend on the eigenvalue. 

To obtain some useful information about the fundamental energy of the radial Hamiltonian p3[) . we need to analyze 
its potential (115[) . The numerator of V2, r(r 3 + a 2 r + 2a 2 M), is negative in some interval to the left of zero, and, 
since Vi(0) = — 4a 2 M < 0, by continuity V\ must also be negative in some neighborhood of r = 0. Thus, there is an 
interval s\(M) < r < S2(M) < where both Vi and V2 are negative. Let ip be a smooth normalized ( — 1) test 

function supported on r G (s%(M), S2(M)). Note that ip has compact support away from r = 0, and that 



(19) 
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FIG. 3: Intersection of the numerically generated curves Ere = 2)(aij)\ <U o=ik and e (fc) for a/M = 1.4. Note the agreement with 
the results in [f|, according to which for a = 1.4 the intersection should occur at k ~ 7.07 x 1.4 = 9, 898. 



with 



82 (M) 

si(Af) 



Vj —= < for j = 1,2. 



Since (^>|Vi|V>) is negative, there is a k c such that, for fc > k c , (tp\ — (d/dr*) 2 \ip) + (ip\V \ip) + k{ip\Vi\ 
and so (■ip\H.\tp) < k 2 (tp\V2\ip) ■ This implies that if —e (k) is the lowest eigenvalue of 7i then 



e (k)< (ip\H\iP) <fc 2 (V>|V 2 |V> <0, iffc>fc c 



From (fTo| and the above equation 

e a (k 



e Q (k) > MVM\k 2 , k>k c 



(20) 
is negative 

(21) 
(22) 



From ([5]), (J6j> and ([22]) it is clear that, for any £ = 2,3, the curves e (k) and Ee,(cuv)\ au> -ik intersect at some ki > 0. 
We have numerically checked this behavior and spotted the predicted intersection point fc^_ 2 ). The results, shown in 
Figure 3, are in perfect agreement with those in [6J. 

This implies that there is a common eigenvalue E = Ei{aui)\ au> —ik t — e (ke) for the radial and angular Teukol- 
sky equations and that {S(9) — S , (^ >m =o,fe <J ) (^), R(r) = At il>i ke \r)}, £ — 2,3,.. are (infinitely many!) solutions 
of the Teukolsky equations, all of which correspond to unstable gravitational perturbations $( s= _2) (t, r, 8) = 

S(g im = k o )(9)A%ipo ke \r)e'xjp(kit/a) with ipo (r) behaving as in l|18[). The radial exponential decay rate guaran- 
tees that any quantity of interest vanishes asymptotically for large r. 

Regarding the behavior of the perturbation at the ring singularity we point out that: (i) Contrary to what happens 
for other naked singularities (e.g., negative mass Schwarzschild solution, see 0,Q)) r = is not a singular point of 
the radial equation, the potential (|15[) being smooth everywhere. In fact, r = is neither a boundary point, nor a 
singular point of the second order radial equation (fT3")) . then imposing conditions at this point, besides requiring that 
the perturbations vanishes at r* — ±oo, would lead to an unnaturally over-determined problem, (ii) It is natural to 
ask, however, what is the effect of the perturbation on the singularity. In particular, one wants to make sure that first 
order corrections to the Riemann tensor invariants do not diverge faster than the zeroth-order term as we approach 
the singularity, otherwise, the use of linearized equations is not justified, since we cannot keep the perturbation small 
in the whole spacetime. As shown in [6|], no algebraic invariant gets corrections under our unstable modes. A deeper 
treatment of the perturbed geometry could possibly be afforded in the approach where the s = — 2 solutions of the 
Teukolsky equation are used as "potentials" to construct metric perturbations [l3l ]. but this is a rather technical 
analysis [19| that certainly exceeds the scope of this report. 

Note that the unstable perturbations are not restricted to the fundamental radial mode, since for large k, the 
potential V in (|13p has a deep minimum that can be approximated by a harmonic oscillator potential of depth of 
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order k 2 , and strength also of order k 2 , so that the level spacing near the ground state of TC is of order k. Thus, the 
negative —e n (k) of the lowest radial eigenvalues grow also quadratically in k for large k, and intersect the angular 
eigenvalue curves for sufficiently large values of k 



IV. UNSTABLE MODES FOR REGION III OF A KERR BLACK HOLE. 

The calculations above can be easily adapted to deal with perturbations in the interior region 
r < r, := M — v 7 M 2 — a 2 of an a < M Kerr black hole. The extreme a = M and sub-extreme a < M cases 
require separate treatments. 



A. Extreme case (a=M) inner region 

The solution of dr* / dr — 1/ \fK in the interior region r < ri = r Q = M is 



( M - r\ 

r* = -lnf ns -J,r<r i , (23) 

with inverse 

r = M(l-e- r *), -oo < r* < oo (24) 
Using the integration factor A 3 / 4 as before, we are led back to ([TBI and (Ti"5|) , with r given in (j2"4"|) . Note that 

^ |4fc 2 exp(2r*) ,r*^oo 

~ |fc 2 exp(-2r*) ,r*^-oo, 1 J 

then the spectrum of the self-adjoint operator TL is again fully discrete and has a lower bound. The eigenfunctions 
behave as 

^K^r «p[- 2fc (*)]( i+ °^) '^ m_ (26) 

[M 1 2 * e S (l + C(M/r)) ,r->-oo 

The argument of instability in the super extreme case goes through in this case without modifications, since the test 
function in (Til?)) is supported in the r < region, and thus can be used again in this case to obtain the bound (|2"Tj) . 
The radial decay (|26[) guarantees that corrections to relevant quantities vanish in these limits. The considerations 
about the perturbation behavior near the ring singularity for the super extremal case apply again to this case, and 
also to the sub-extreme case below. 



B. Sub-extreme case (a < M) inner region 

The sub-extreme case introduces some subtleties: requiring dr* /dr = 1/vA for r < rj gives 



* . / M -r- vV 2 -2Mr + a 2 \ n+r -2r- 2 v /(r - r)(n -r)\ 
r =ln [ M ) =ln [ — a J' r<n > (2?) 

and the radial equation reduces to P^|) . (fTS)) with r the inverse of J2"7|) . 
Since 

.* ^ „* i„ i r ° — r « 



oo < r* < r* := In , (28) 
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(I13p is, in this case, a Schrodinger like equation on a half axis, with potential diverging as 

V ~ [k(M 2 - a 2 )/(2aM)} 2 exp(-2r*) 

for r* — > — oo and, for r* — > as 



V 



1 4Af 2 (M - VA/ 2 - a 2 ) 



VM 2 



M 2 



a M (M - VM 2 - a 2 ) 



1 



2 _ 1 
4 



^7** 7 1 * ^2 



(29) 



> 0. The subtleties mentioned above are the half infinite domain for r* and the behaviour (f2"9"]l of the potential 
at the r* boundary. These singular potentials are discussed in detail in (2fj| . The eigenfunctions of 7i which are square 
integrable for r* near — oo behave as in (fT8| in this limit and as 

rf> ~ a [(r* - r*)3+" + ...] + 6 [(r* - r *)*-" + ...] , (30) 

near the horizon, the E eigenvalue appearing only in sub-leading terms. Thus, if v > 1, for generic E these are not 
square integrable near zero, unless 6 = 0, and this is precisely the condition that selects a discrete set of possible E 
values as the spectrum of H, and that defines de space of functions where Tt is self-adjoint. On the other hand, if 
v < 1, for any E the eigenfunction behaving properly at minus infinity will be square integrable for r* € (— oo,0), 
and a choice of boundary condition needs to be imposed to define a set of allowed perturbations D p h ys , in order that 
7i be a self-adjoint operator on D p h ys (and thus have a complete set of eigenfunctions) [2l| . This situation, where 
a choice of self-adjoint extension needs to be done, is also found in the study of perturbations around the negative 
mass Schwarzschild space-time Q, and perturbations in the inner region of a Reissner-Nordstrom space-time fl9j [. 
Roughly, Dphys is constructed by adding to the set of functions of compact support, those decaying as (f3"0")) with 
a fixed quotient 5 := d/a 6 R, possibly infinite [2l[. However, regardless of our choice for q, functions of compact 
support away from this singular point will belong to D p h ys , in particular, the test function used in (|19p -(|20 p . This 
implies that the bound (|21|) will hold for any choice of boundary condition (q value)! In other words, region III 
of a sub-extreme Kerr black hole is unstable regardless the boundary condition imposed to perturbations at the 
inner horizon. A natural choice would be that giving the fastest decay, i.e., q = 0. In this case the gravitational 
perturbations behave as $ ~ A 3 / 4 (r* — r*)^ +v ~ (r* — r*) 2+u ~ (r, — r) 1+ly / 2 . One should keep in mind that 
v = v(ki) is positive and grows for large I. 



The instability found in region III of Kerr space-time, and the fact that the Reissner-Norsdtrom charged black 
hole also has a two horizon structure with an inner static region r < r i: triggers the question of whether of not this 
inner static region is stable. Preliminary work indicates that this region is unstable [191 ]. 

As a final comment, we note that our results do not contradict the established stability of the outer stationary, 
region I of Kerr black holes. For a < M and r > r Q we can define r* as in (fTTj) . and get again the Schrodinger form 
(|13[) . but now the domain is restricted to r D < r, where 



V > 



15 

T 



if a < M and < k. 



(31) 



This implies that e {k) < 15/4 for all k > and thus cannot intersect E g{auj)\ aU i=ik for any t. To get the bound ([31 
is easy in the extreme case, for which (|15|) reduces to (replace r — Mr' , then drop the primes) 



Vr. 



r(r + l)(r 2 - r + 2) 





'r 2 - 2r- 1" 




k 2 -4 








r- 1 





(r- l) 2 

We search for critical points at fixed k by solving dV/dr = 0. This gives 

2(r-l)(r 2 -2r + 3) 
~ ( r 2 + l)(r 2 -2r-l)' 



(32) 



(33) 



which has a unique r > 1 = r$ = r Q solution for every positive k, with r moving from infinity down to 1 + v2 as k 
goes from zero to infinity. Since d 2 V/dr 2 > at these points, these are (global) minima of V for fixed k. Moving 
along this curve from r = 1 + v2 to infinity, V decreases monotonically to —15/4. In the sub-extreme case a < M 
the equation dV/dr = is quadratic in k, and yields two curves like that in ([33]) . However, one of these curves gives 
k < for all r > r ext , and should be discarded. The other one corresponds to minima of V for fixed positive k, with 
V again decreasing monotonically to —15/4 along it. 
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